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The intrinsic spin-orbit coupling (ISOC) in bilayer graphene is investigated. We find that the 
largest ISOC between tt electrons origins from the following hopping processes: a tt electron hops 
to G orbits of the other layer and further to tt orbits with the opposite spin through the intra- 
atomic ISOC. The magnitude of this ISOC is about OAGmeV, 100 times larger than that of the 
monolayer graphene. The Hamiltonians including this ISOC in both momentum and real spaces 
are derived. Due to this ISOC, the low-energy states around the Dirac point exhibit a special spin 
polarization, in which the electron spins are oppositely polarized in the upper and lower layers. This 
spin polarization state is robust, protected by the time-reversal symmetry. In addition, we provide 
a hybrid system to select a certain spin polarization state by an electric manipulation. 



Recently, ultrathin graphitic devices including mono- 
layer and bilayer graphene have been fabricated 
experimentally[Tl [2]. The experimental developments 
greatly stimulate theoretical research on this subject [3] ■ 
Monolayer graphene is a gapless semiconductor with 
unique k-linear and massless Dirac-like spectrum, while 
bilayer graphene presents parabolic and massive spec- 
trum. This distinction in dispersion leads to many physi- 
cal diversities. To name a few, the Hall plateaus in mono- 
layer graphene appear at half-integer values while that 
in bilayer graphene appear at integer values without the 
zero- level plateau [2 [21 HHB]. The monolayer graphene 
exhibits the weak anti-localization behavior due to the 
TT Berry's phase, however, in bilayer graphene the 2n 
Berry's phase leads to weak localization |7J. The differ- 
ence between monolayer and bilayer graphene has led to 
great effort in research. 

The spin-orbit coupling (SOC) has attracted great at- 
tention in the last decade as it plays a very important 
role in the field of spintronics. SOC can couple the elec- 
tron spin to its orbital motion and vice versa, thereby 
giving a useful handle for manipulating the electron spin 
by external electric fields. Furthermore, the SOC may 
also generate some exciting phenomena, e.g. the spin 
Hall effect and the topological insulator. [8| The SOC ori- 
gins from the relativistic effect, so it is usually very weak. 
For example, in a monolayer graphene the theoretical re- 
search both in tight-binding model and the first principle 
calculation show that the strength of SOC there is only 
of the order 10~^eF.[9] It will be very interesting to find 
a large SOC or to enhance the SOC in a graphene system 
by some means. 

In a carbon atom, the SOC AL ■ a (with the angular 
momentum L and spin a) can couple the 2pz orbit (or tt 
orbit in the graphene) to the opposite-spin 2px and 2py 
orbits (or a orbits in the graphene). The magnitude A 
can be estimated to be of order 5meV. In monolayer 
graphene, a orbits hybridized in a sp^ configuration to 
form the honeycomb lattice, while, the remaining 2pz or- 
bit, which is perpendicular to the graphene plane, forms 
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FIG. 1: (color online) Schematic diagram for hopping pro- 
cesses in bilayer graphene, which leads the largest ISOC be- 
tween TT electrons. The black (red) point stands for A (B) 
sublattice, respectively. 



the TT orbits and dominates the low energy physical prop- 
erties of graphene. Because the tt orbits are perpendicu- 
lar to the sp^ orbits, the hopping process from tt orbits 
to all a orbits are suppressed. Therefore the intrinsic 
SOC (ISOC) in the low-energy tt electrons must involve 
three hopping events: the SOC between the spin-up tt 
orbit and the spin-down a orbit, the hopping coupling 
from the a orbit to its nearest-neighbor (N.N.) a orbit 
with the same spin, and the SOC from the spin-down 
a orbit to the spin-up tt orbit in the N.N. atom. This 
ISOC is on the order and its strength is very small, 
about 10~^ey,[n]. There are several ways proposed to 
enlarge the SOC strength. pUl [TT] e.g. by mixture of tt 
electrons with a electrons on the N.N. sites by introduc- 
ing impurities, but this approach will cause the system 
to be unclean. 

In this Letter, we predict a large ISOC in a clean bi- 
layer graphene. We consider bilayer graphene which con- 
sists of two coupled honeycomb lattice with sublattice 
Al, Bl (A2, B2) on the top (bottom) layer, respectively. 
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arranged by Bernal stacking as illustrated in Fig. 1. Ev- 
ery Al site lies directly on the top of A2 site, while Bl 
and B2 do not. We notice the bond of B1-B2 is not per- 
pendicular to sp^ orbits, so the coupling between the tt 
and a orbits of the inter-layer B atoms can exist with the 
strength on the order of Q.leV . Then as illustrated in Fig. 
1, the TT electron (red circle) at site Bl can hop to a or- 
bits (green circle) at site B2 and further hops to the B2's 
TT orbit with the opposite spin through the intra-atom 
ISOC. This hopping process contributes to the largest 
ISOC between the tt electrons in bilayer graphene, with 
its strength being on the order of A (not A^), so this 
ISOC should be much large than ISOC in a monolayer 
graphene. 

In order to derive the ISOC between tt electrons in bi- 
layer graphene, we need to consider one s orbit and three 
p orbits for each site, together with two spin indexes, two 
sub-lattice indexes and two layer indexes, totally there 
are 32 orbits in every cell. Then, the tight-binding Hamil- 
tonian reads: 



H 



(1) 



where {^l) is the annihilation (creation) operators of 
an electron in the cell j (i) with its dimension 1 x 32 (32 x 
1). iJjri Hg-, and J7 are the Hamiltonian of tt electrons, 
a electrons, and the coupling between tt and a electrons. 
We adapt the parameters of SWM model, [HHH] and the 
TT electron Hamiltonian H^^ includes the intra-layer N.N. 
hopping coupling 70, the interlayer hopping coupling 71 
of A1-A2 and that 73 of B1-B2 (see Fig.l). The value 
of 70 = 3.12eF,7i = 0.38eF,73 = 0.29eV is adopted 
from Ref. |12l[Tl] . The on-site energy of tt electrons is 
set at zero as an energy zero-point. For the cr electrons, 
we choose the atom orbits 2px, 2py, and 2s as the base 
vectors. The on-site energies of Hamiltonian iJ^- are 
for 2px and 2py orbits and esp {cgp = —8.87eV) for 2s 
orbit. There is no coupling between a electrons in same 
atoms because of the orthogonal base vectors. But the 
hopping couplings between intra-layer N.N. A-B atoms 
are very strong. We first consider a special case that the 
bond of A-B atoms is parallel to x axis, thus the hopping 
Hamiltonian H^h reads: 



ah,AB 
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where ^aA/B = {a2px,A/B,a2py,A/B,a2s,A/BV is the an- 
nihilation operators of electron in the 2px, 2py, and 
2s orbits of the A/B atom. In Eq.(2) the B atom 
must be the N.N. of A. The hopping elements Vppcr — 
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3.12ey, Vs, 



5.58ey, and 



is at an angle 9 to the x axis, by taking a rotation in 
xy plane, the hopping Hamiltonian H^h is obtained as: 
Hah = ^AR~HO)VhopR{0)^lB + h.c, with 



m = 
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The hopping between a orbits of inter-layer is neglected 
because that the a electron is mainly within each layer. 

The TT-a hopping Hamiltonian J' is from both ISOC at 
each atom and the inter-layer ir-a coupling. The intra- 
layer TT-a coupling for the N.N. A-B atoms is zero as 
mentioned above. The inter-layer tt-<t coupling in the 
aligning Al and A2 atoms is also zero since the Al's 
TT orbit points to the center of the A2 atom. But the 
inter-layer TT-a coupling between the Bl and B2 atoms 
exists because that the bond B1-B2 is not perpendicular 
to the 2px and 2py orbits. Such kind of TT-a coupling is 
a crucial characteristic in bilayer graphene and does not 
exist in a monolayer graphene. The Hamiltonian of this 
TT-a coupling is: 
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+ h.C., (4) 



where 73 = 73 * a/d w O.leV with the lattice constant a 
and the distance d between the layers. is the angle of 
the projection of the B1-B2 bond in xy plane and the x 
axis. For every Bl atom there are three N.N. B2 atoms, 
and their TT-a couplings are similar except for the different 
(j). Similarly, the B2's tt orbit also has the coupling with 
the a orbits of the Bl atoms. Finally, in the 2s and 2p 
base vectors, the SOC AL ■ a at each atom site can be 
easily rewritten as: 



Jsoc 



— fi ' 

cy pza 



^pya 



h.c. 
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where a =t, i is the spin index, a —\, (f) for a =t (i) j 
and i is the imaginary unit. Here the SOC is the spin- 
flip hopping, but all other hoppings mentioned above are 
between the same spin orbits. 

Because that the infinity bilayer graphene is a periodic 
system, we can rewrite the Hamiltonian in Eq.(l) in the 
momentum space: 



^ppa 





Vsp 









^pp-K 





ii^B + h.c. 
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Vsp 





Vss 







{dk,dky/4TT^)¥^ 

{dk,dky/4TT^)i'lH{k)i'k, 



(6) 



where the momentum k = (k^, ky) is in the first Brillouin 
Zone. 

Next we focus on how to derive the low energy effec- 
tive Hamiltonian. First the high-energy a electrons need 
to be eliminated. For this purpose, we perform a canon- 



Vss = 6.77ey.[TS] Second, if the bond of A-B atoms ical transformation: 'H(k) 



'77(k)e'5,P [ig where 
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FIG. 2: (Color online) (a) and (b) are the dispersion rela- 
tions for infinite bilayer graphene (set ky = 0) and zigzag 
bilayer graphene nanoribbon with width W = 120\/3a. (c) is 
the schematic diagram for a hybrid device consisting of two 
large bilayer graphene connected by a narrow zigzag bilayer 
graphene ribbon. The green (blue) arrows in (c) represent the 
spin polarized direction for the green (blue) state in (a). 



S = 



is anti-hermitian matrix. After the 





-A/t 

canonical transformation, we demand that the a and tt 
electrons are decoupled in the Hamihonian H(k). By 
using the condition of the decoupling and keeping the 
matrix elements in up to order J^^, the matrix M 
can be obtained: 



M 



— "^J H„ 



(7) 



While around Dirac points for the low-energy electron, it 
is safe to take the approximation M « —H~^J', because 
of large energy gap of several eV between a and tt elec- 
trons. To substitute this M into the above equation of 
the canonical transformation, the effective Hamiltonian 
of TT electrons with ISOC reads: 



(8) 



The magnitude of the ISOC term is roughly estimated as 
^Ts/^o-TT, in which ecnr is the energy difference between 
a and tt bonds at Dirac points. 

In bilayer graphene, the energy bands of tt electron 
are fourfold except for the spin freedom of degree. Two 
energy bands touch at Dirac points dominate the low- 
energy properties, while the other two bands have an 
energy gap of about 271 due to the strong interlayer cou- 
pling between tt electrons on site Al and A2. Because 
that the gap 271 is about 0.7eV, we need to carry out 
the canonical transformation procedure again to elimi- 
nate the high energy tt bands for sites Ai and ^2- Fi- 
nally, we get the low-energy effective Hamiltonians Hk 



and Hk' around the Dirac points K and K': 
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where Z = { ^^^^ -I- 373afc_}, vp ~ 37oa/2 is the Fermi 



- iky J k— — kx ' 



velocity of tt electrons, and fc-|_ = k^ 
In the Hamiltonians Hj^ or H^', the momentum k is 
measured from K or K' , respectively, and the basis or- 
bits are (p^^Bl, pz^Bl, pz^B2, pziB2). The first term 
in Eqs.(9) or (10) comes from the well known low-energy 
Hamiltonian of tt electrons on sites Bi and B2 of bi- 
layer graphene. [51 [T^ The second term is the ISOC be- 
tween the TT electrons, which is a central result of this 
work. This ISOC origins from the interaction of both 
the intra-atom ISOC and interlayer hopping that cou- 
ples the TT-a orbits as shown in Fig.l. Its strength m is 
about Ay^/ea-K, which is on the order of A. From our 
calculation, the strength m = 0A6meV and it is about 
100 times larger than the ISOC in a monolayer graphene. 

Eqs (9) and (10) lead to energy bands E{k) — 
±1/2(|to| ± -f 4|Z|2). Fig.2a depicts the energy 

dispersion at fixed ky = 0. Due to the ISOC, the spin 
degeneracy is removed. Here two bands have an energy 
gap 2m and the two low-energy bands still touch at Dirac 
points. Around the Dirac point K, the low-energy bands 
are spin polarized. The spins in the upper layer are 
mainly upward while they are downward in the lower 
layer, as shown in Fig. 2c. Moreover, this spin polariza- 
tion is almost independent to the momentum |k|. To 
compare with the quantum spin Hall effect, in which the 
spin polarization is opposite on the opposite edges, now 
the opposite spin polarizations exhibit on the opposite 
surfaces. Near the Dirac point K' , the spin polarization 
in the low-energy excitation mode is in contrary to that of 
K point since they are the time-reversal Krammer pair. 
Given a time-reversal invariant impurity, the inter-valley 
scattering processes between these two Krammer's pair 
are suppressed. jlTj il8j Thus, the spin relaxation time is 
long. 

After obtaining the ISOC term in k-space, its tight- 
binding form in real space can be derived straightfor- 
wardly. Here we keep the tt orbits of the Al and A2 
atoms, and the tight-binding Hamiltonian is; 



(i+<5,)c 
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where aia(al^), Cia(q^), and dia{dlj denotes 

the annihilation(creation) operator of the tt electron at 
the Al, Bl, A2, and B2 atoms with spin a in the cell 
i. "(?j=i,2,3 are the three N.N. vectors. In equation (11), 
the first, second, and third terms express the intra-layer 
and inter-layer 7r-orbit coupling without spin flip. The 
last term leads to the spin flip hopping between the Bl's 
and B2's tt orbits due to the ISOC, with tr — t% — 

^exp{iK - Si). We note that this Hamiltonian maintains 
C3 symmetry as well as the time reversal symmetry. By 
using this tight-binding model, electronic properties of 
finite bilayer graphene can be calculated. 

Until now, we exhibited that the low-energy modes 
in bilayer graphene are spin polarized due to the ISOC. 
But in the equilibrium both modes around K and K' 
are occupied and the total spins are unpolarized. For 
the purpose of selecting a single excitation mode, we 
need to choose a valley. We design a hybrid system, the 
so called 'valley filter' [19] consisting of two wide bilayer 
graphenes connected by a narrow zigzag bilayer graphene 
nanoribbon (see the sketch map in Fig. 2c), to select a sin- 
gle excitation mode. Earlier works [TOl [2U] on monolayer 
graphene ribbon with zigzag edges have shown that there 
exists a propagating chiral edge state, and the 'chiral- 
ity' means that the sign of group velocity determines the 
valley. The situation is the same for the zigzag bilayer 
graphene ribbon, as demonstrated from its dispersion re- 
lation in Fig. 2b. Let us set the Fermi energy as shown in 
Fig. 2a and 2b and apply a positive voltage to the device. 
There will be two kinds of excitation modes (blue and 
green arrows in Fig. 2c) in the left region. However, the 
edge states of the central zigzag ribbon can only support 
the one around K point, thus only one excitation mode 
(green arrows in Fig. 2c) transport to the right region 
after the selection, which is spin polarized. Protected 
by time reversal symmetry and the large separation be- 
tween K and if, the spin polarization in the right region 
is difficult to be scattered and can last for a long time. 
Meanwhile, if we change the direction of current flow, 
the other kind of excitation mode can present in the left 
region. 

In summary, by considering the hopping coupling be- 
tween interlayer ir-a orbits and the ISOC in carbon 
atom, we find a large ISOC between tt electrons in bi- 
layer graphene, with its strength being around 0.46meT^, 
about 100 times larger than that in the monolayer 
graphene. The Hamiltonians of this ISOC in both mo- 
mentum and real space are derived. Due to the ISOC, 
the low-energy electrons around the Dirac point exhibit 
a special spin polarization, in which the electron spins 
in upper and lower layers are polarized in opposite direc- 



tions. In addition, we propose a hybrid system to select a 
certain spin polarizing mode by an electric manipulation. 
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